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Abstract 

High precision analytical approximation is proposed for variance-covariance based risk allocation 
in a portfolio of risky assets. A general case of a single-period multi-factor Merton-type model with 
stochastic recovery is considered. The accuracy of the approximation as well as its speed are compared 
to and shown to be superior to those of Monte Carlo simulation. 



1 Introduction 



Merton-type model s, also referred to as structural models, such as PortfolioManager (jKealhoferl . 120011 ) 
and CreditMetrics (jGupton et all I1997T ) . have become a standard choice for financial institutions' credit 



risk economic capital frameworks. In these models default correlations between different borrowers are 
modeled using a set of common systematic risk factors associated with the state of economy. Computa- 
tionally heavy Monte Carlo simulations are usually used for calculations of portfolio-wide risk measures 
as well as risk allocation to sub-portfolios and/or individual exposures. However, simulation-based risk 
allocation on exposure level suffers from Monte Carlo noise and is especially demanding in terms of 
computer power/time. 

Advanced Monte Carlo simulation techniques have been developed in order t o improve the stability 
problem associated with risk decomposition down to individual facilities (see, e.g.. lTasche[|2009l ). Despite 
the improvements, the problem exists due to the inherent stochastic nature of Monte Carlo simulations. 

Practical demand for stable and fast capital allocation routine in credit portfolios led to development 
of analytical techniques. Analytical allocation techniques are also preferable for the purposes of portfolio 
optimi zation and risk adjusted pricing. Granularity adjustment ( Gordv . 20041 ) and multi-factor adjust- 
mgnt (|Pvkhtml . l2004 ) are well known extensions of asymptotic single risk factor framework of I Vasicekl 
(|2002h . 

Unfortunately, most of the research on risk allocation techniques focuses on 'advanced' risk measures 
like VaR and ES (expected shortfall), leaving va riance-covariance based a llocation approach aside. Despite 



the shortcomings of this approach (see, e.g., Kalkbrener et al. . l2004h . this (old-fashioned) allocation 



method still remains the allocation method of choice for many financial institutions. Variance-covariance 
ba sed risk contr ibution is considered to be intuitive, simple and easy to compute risk measure (according 
to lBCBSl . [2009h . Yet, no efficient analytical solution has been reported so far. A brute force approach 
consists of calculations of all pairwise correlations in the portfolio (the reader is referred to Section [2~2l 
for the details). The amount of such calculations is quadratic in the number of loans in the portfolio. 
This quadratic complexity of the calculations make such an approach impractical for big portfolios. 

In this article a variance-covariance based analytical risk allocation technique is proposed. The pro- 
posed approach is applicable to fully featured Gaussian multi- factor Merton-type models, is suitable for 
virtually any portfolio size and composition and is remarkably accurate and fast. The main advantage of 
the proposed technique is that the underlying algorithm is of linear complexity in portfolio size. 
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It should be pointed out that the author is not an advocate for a particular risk measure. Financial 
institutions have the responsibility for choosing the risk allocation approach and should clearly under- 
standing the associated shortcomings. The outcome of any risk allocation technique should be a subject 
of regular expert reviews and sanity checks. 

The technique presented here might be of interest to practitioners building in-house economic capital 
models. In contrast to most publications focusing on default-only case, the proposed approach is suitable 
for a very general class of credit portfolio models. 

The article is organized as follows. First, a general case of a portfolio of risky assets, whose pairwise 
correlations are modeled by means of a set of normally distributed common factors, is considered. A 
routine to allocate portfolio's risk (associated with its standard deviation) to underlying exposures is 
described. Next, it is demonstrated how to apply the proposed allocation technique to PortfolioManager- 
like credit portfolio model, i.e. single-period, multi-factor Merton-type model with stochastic recovery 
and risk neutral valuation at horizon (which accounts for credit migrations). Finally, accuracy and speed 
of the proposed technique are compared to those of Monte Carlo simulation. 



2 Theory 

In this section a formal theoretical framework for variance-covariance based risk allocation is described. 
First, it is shown how the Euler allocation principle is applied to portfolio's standard deviation and a link 
with risk allocation in credit portfolio is briefly described. Next, Gaussian multi- factor model is considered 
and a series expansion for pairwise covariance is proposed. Finally, an algorithm for calculation of the 
risk contributions is described. 



2.1 Standard deviation and risk allocation 



Standard deviation is a subadditive risk measure allowing Euler allocation (see, e.g., TascheL 20081 ) of 



portfolio's risk to sub-portfolios and/or individual exposures. Consider a portfolio consisting of risky 
assets with standard deviations {ci} and pairwise value correlations {pij}. The standard deviation of 
the portfolio o~ p , 

can be written as a sum of Euler risk measures (risk contributions) o~f as 

i i 

The risk contributions {of } can readily be calculated by differentiation and can be expressed in terms of 
covariance (or correlation pi p ) between individual exposures and the portfolio 

„ cov(z,p) 1 v , . , . 

= °iPip = = — V c «v(i, j). (2.3) 

CTp CTp 

The variance-covariance approach to capital allocation in credit portfolios utilizes the above decom- 
position to distribute the economic capital of the portfolio between individual exposures proportionally 
to the risk weights (|2.3|) as follows 

o c 

(capital charge)i = — x (total economic capital). (2-4) 
cr p 

2.2 Variance-covariance structure of Gaussian multi-factor framework 

Let us start by considering a portfolio of risky instruments {vi}. The value u, of each instrument is 
assumed to be a function of normally distributed random variable e^. Correlations between variables {e^} 
are modeled through a set of normally distributed independent variables {r]k} referred to as common 
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factors. Each variable q is split in a sum of instrument specific part, which depends on a Gaussian 
variable £i, and systematic part, which depends on the common factors, as follows 



ti = n Y,((3i)kVk + ^l-rUi- (2-5) 

A; 

The independently distributed random variables {{£i}, {?7fc}}0 are assumed to have zero mean and unit 
variance. Instrument specific constants and {({3i)k} determine dependency of q on the common factors. 
The so-called factor loadings {{Pi)k} are subject to normalization condition 

£(A) 2 fc = I- (2.6) 

k 

Using the representation (|2.5|) . one can readily calculate a covariance p^ between any pair of distinct 
variables {e^} 

Pij = (eitj) = UTj ^{Pi)k{(i 3 )k = nrj^iPj, i ^ j. (2.7) 

k 

The value covariance between two distinct instruments (viVj) can be written as 

(viVj) = Jvi(€ i )v j (€ j )n 2 (€i,e j ,pi j )deide : j - Jvi(€i)n(€i)dei j " Vj(ej)n(ej)dej, i^j, (2.8) 

where n and n 2 are normal and bivariate normal density functions defined as 

n(x) = -±=e-* 2 '\ n 2 (x, y, p) = - 1 2 exp (- *- ' ~? pX% . (2.9) 
V2tt 2tt^1 - p 2 V 2(1 - p 2 ) J 

In theory, the above expression and (|2.3p are sufficient to calculate the standard deviation based 
risk contributions {erf}. In practice, however, one is facing a problem of computing N(N — l)/2 two- 
dimensional integrals in (|2.8[) for a portfolio of N instruments. Such a brute force approach becomes 
inappropriate for large portfolios, i.e. N = 10 4 or higher. 

2.3 Series expansion for covariance 

Aiming for a linear (in N) complexity algorithm for computing the risk contributions ()2.3[) , let us reduce 
the d imensionality of the integral in (|2.8p . Applying Mehler's formula (for the proof see, e.g., Foatal . 



V*TT f Vff < \P 1 ( 2pxy-p 2 {x 2 +y 2 ) \ 

> He„ (x He„ (y)— = exp — - (2.10 

n\ Jl-p 2 V 2(1 -p 2 ) J 



19781) 

2 / 2 

— mm I 

- V^J 2 

to the bivariate normal density function, the following expression is obtained for the pairwise value 
covariance 

oo n n 



p n f f 

= z2^f fi(ei) He n(e0 n ( e de i / u i( e i) He n( £ j) n ( e i)de^ (2.11) 

n=l H ' J J 

where He n (x) = (— l) n e x2 / 2 (d/dx) n e~ x2 / 2 are Hermite pol ynomials (for definition and properties of 
Hermitc polinomials see, e.g.. lAbramowitz and Stegu l ll972l ). The last expression can be rewritten as 

(ViVj) = E PiA n)v l"\ vin) = ^J v(e)Re n (e)n(e)de. (2.12) 



The above series conver ges at least as fast as a geomet ric series with quotient pij. To see this, one 
can use inequality (see, e.g., Abramowitz and Stegunl . 1972h 

Ke n {x) < kVri.e x2/4 , few 1.086435, (2.13) 



1 Assuming all these variables to be independently distributed is equivalent to an assumption that each borrower in the 
portfolio is represented by one facility. This assumption will be relaxed in Scction l3.2l 
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which leads to 

\ V W\ < 4r / k(^)l|He„(e)|n(e)de < -±= f \v{e)\e~^ ^de. (2.14) 
Vfi! J V27T j 

The last integral in the above inequality is finite for any reasonable value function v{e). 
2.4 Linear complexity algorithm 

Using the series expansion (|2.12p and expression (|2.7[) . a linear complexity algorithm can be constructed 
for calculation of risk contributions (|2.3|) as follows. 

(n) 

First, the parameters v\ , n = 1 . . . n maa: are evaluated for each instrument i = 1 . . . N. The number 
of terms n max used for calculations determines an accuracy of the algorithm. Second, assuming Nf 
common factors, the following portfolio specific parameters are calculated 

N 

P £Lk n = E r ? w i n) (A)fci • ■ ■ (AK, % = 1 • • ■ (2.15) 
Third, the risk weights f|2.3|) are calculated using the following expression 

*x = <«?) + E < v ^> . ( 2 - 16 ) 

which after substituting (|2.12p . (|2.7p and (|2.15p becomes 

a p at = a*+Y: r?«W £ ^ . . . (ftJ^pW^ - £ ( r? „«) . (2.17) 

n.=l fei...fc n n=l 

Finally, the standard deviation of the portfolio and the risk charges per instrument are calculated as 



^ = JE<V^ °t= * 1 c - (2-18) 

Under certain circumstances the algorithm may be further simplified. Consider, for example, a single 
factor model. Expressions (|2.15|) and (|2.17p then become 

"max "max 2 N 

= + E r?4 n) p^ - E { r >i n) ) ■ = E tf«< (w) - ( 2 - 19 ) 

n — 1 71—1 i—1 

Another simplification takes place if we consider a default-only credit portfolio, i.e. 

In this case the parameters can be computed analytically 

e" d2 /2 

«<»> = Z^_He n _x(d) (2.21) 

3 Practice 

In order to assess the accuracy of the proposed risk allocation routine, PortfolioManager-like credit port- 
folio model is considered. In this section basic description of the model is given. Some aspects concerning 
application of the proposed risk allocation framework to the credit portfolio model are discussed. The 
analytically calculated risk measures are compared to those obtained as a result of Monte Carlo simula- 
tion. It is finally concluded that the accuracy of the analytical approximation is higher than the accuracy 
of Monte Carlo simulation with 10 8 scenarios. 
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3.1 Credit portfolio model 

Consider Gaussian multi-factor Merton-type credit portfolio model with stochastic recovery rate. The 
credit portfolio consists of loans whose values {vi} are functions of the corresponding borrowers' asset 
returns {ti}. The asset returns are assumed to be normally distributed with the correlation structure as 
described in Section 

A borrower is considered to be in default at some point in time if its asset return falls below the 
borrower-specific threshold, i.e. when e < —d. The quantity d is called distance to default and is related 
to the borrower's probability of default p through the standard normal cumulative distribution function 
$ as 

f-d e -e*/2 

>-«-*>- Lust*- (3 ' 1) 



The portfolio is observed at some time in the future th called horizon. Assume that both probability 
of default till horizon p and (cumulative) probability of default till maturity p rn are known for each loan 
in the portfolio. In case a loan matures before the horizon, its probability of default at horizon is assumed 
to be equal to its probability of default at maturity. 

In case of default, the value is defined by the recovery rate (1 — and loan's risk-free value VQe~ r ^ tm ~ thS> 

as 

v(e) = (1 - l)voe- r ^-^ , e<-d, (3.2) 

where r is a risk-free rate and t m is loan's maturity. In case of no default and maturity before or at 
horizon, the risk-free value of the loan is recovered 

v(e) = voe- r ^- th \ t m <t h , e>-d. (3.3) 

Risk neutral valuation, as described in I Vasicek (l2002h . is applied in no default case to the loans maturing 
after horizon 



where 

b = <f- 1 ( Pm ) + Xr t -^ w ^. (3.5) 



The A in the above expression stands for the so-called market price of risk and r is the borrower-specific 
parameter introduced in (|2.5p . The above valuation takes into account change of loan's market value due 
to credit migration. 

Uncertainty in recovery rat es in case of default is modeled by means of Beta distribution (see, e.g., 



Abramowitz and Stegunl . 119721 ) whose probability density function is 

f( X> a,b) = — L-x-^l-x) 6 - 1 , (3.6) 
B(a, b) 

where a and b are parameters characterizing the distribution. The above parameters are chosen to math 
the mean I and the variance a 1 of the loss distribution 

a+b a+b+1 k K ' 

where k is a parameter defining the shape of the distribution. Losses in case of default are assumed to 
be distributed independently from asset returns. 

3.2 Applying risk allocation routine 

In order to adopt the risk allocation procedure proposed in Section 12.41 some modifications should be 
made to (f2~T6"|) and (f2~l7]) . 
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For the value function v(e) of a loan at horizon which is defined by (|3.2j) - (|3.4|) . the standalone variance 
a 2 has to be calculated as 



a 2 = Jv 2 (e)n(e)de- (J v{e)n(e)de\ + (yoe-^-^Y , (3.8) 

where the last term takes into account loss uncertainty in case of default according to (|3.7[) . Since 
distribution of losses in case of default is assumed to be independent, the other terms in (|2.17p are not 
affected by stochastic nature of recovery rate. 

Expressions (|2. 16|) and (|2.17| were derived assuming that all the borrower-specific parts of asset returns 
in (|2.5p ) are independent. In practice, however, some borrowers are represented by multiple loans in 
the portfolio. Such loans have perfectly correlated asset returns, which makes the series expansion (|2.12D 
invalid. It is still possible to modify the expressions (|2.16p and (|2.17p . preserving the linear complexity 
of the risk allocation algorithm. 

Consider a borrower a with multiple loans {vi a }. Expression (|2. 16|) has to be modified as follows 

= E ( U ^Vja) + E ( V ™ V jb) • (3-9) 

Introducing V a (e a ) = Vi a (e a ), net value of the loans of the same borrower, the first term in the above 
expression can be written as (vi a V a ). Assuming the recovery distributions of the loans of the same 
borrower perfectly correlated, the following modified version of (|2.17| can be derived 

a v at = (v ia V a ) + £ rrf £ (#)fa ■ ■ ■ (A)*»^?..* B - £ ^iaV^, (3.10) 

n=l ki...k n n=l 

where the first term is defined as 

(vV) = J v(e)V(e)n(e)de - J v(e)n(e)de J V{e)n{e)de + p 1 ^ ~ l K V e- 2r ^- th l (3.11) 

3.3 Numerical results 

To assess performance of the proposed approximation, analytically calculated risk contributions have 
been compared to those calculated by means of Monte Carlo simulation. The comparison has been made 
using the settings described in Section l3~H In particular, the same valuation function (|3.2|) - (|3.4p has 
been used in both cases. Analytical risk contributions were computed numerically using the algorithm 
described in Section [2~4l modified by (|3.10p and (|3.1ip . Monte Carlo risk contributions were estimated 
using the definition (|2.3p according to 

j l N mc _ _ N mc 

at = — ■ irr- ^(vtk - v l )(v pk - v p ), <7 p = ^2(v pk -v p ) 2 (3.12) 
<J P iv mc fc=i fc=l 



and 



N mc „ N. 



— ^2v ik , v p = -^—^2vp k7 (3.13) 



k— 1 K— 1 



where N rnc is a number of scenarios and Vik and v p k are loan and portfolio value realizations in the fcth 
scenario. 

A "real life" portfolio was used for numerical tests. The portfolio consisted of 8036 loans to 4378 
borrowers. Horizon th was chosen to be 1 year, 4% flat risk-free interest rate was used and market price 
of risk A was assumed to be 0.4. The recovery uncertainty parameter k in (|3.7p was set to be the same for 
all loans and equal to 4. A wide range of loan risk drivers was covered. Namely, probabilities of default 
at horizon were in the range from 10~ 5 to 0.4, maturities varied between 1 month and 30 years, losses 
given default I were spread between 0.1 and 0.99. The borrowers represented a wide variety of geographic 
regions and industries. The borrower-specific coefficients {r 2 }, determining dependency on the common 
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Figure 1: Relative differences between Monte Carlo and analytical esti- 
mates of the risk contributions a°. 



factors {?7fc}, covered the range from 0.07 to 0.65. The common factor model used in the experiments 
contained 120 factors {rj^}. 

The accuracy of the Monte Carlo based estimate of a c depends on the number N of scenarios in a 
simulation. The standard error of the Monte Carlo estimates is proportional to . Several Monte Carlo 
simulations were performed with different numbers of scenarios. Starting at 10 5 , the number of scenarios 
was gradually increased up to 10 8 . 

The accuracy of the analytically estimated er c , on the other hand, is limited by the number of terms 
n-max in the series expansion (|2.12[) that are taken into account. The amount of common factors Nf = 120 
in the model imposed a practical limitation on the cutoff point n max in (|3.10|) . since for the nth order 
terms in the series expansion (|2.12j) one has to keep track of {Nf) n portfolio parameters k . Two 
values of the cutoff parameter n max = 2 and n max = 3 have been considered in the experiments. 

Deviations between Monte Carlo and analytical estimates of risk contributions, gr " 1 ^7°' a " , have been 
studied for different combinations of N and n max . Figure [T] displays the relative differences between the 
Monte Carlo and analytical estimates of the risk contributions for each loan in the portfolio together 
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Figure 2: Comparison between the accuracies of Monte Carlo simulations and 
analytical approximation. 



with the standard error of this deviation. One can see that when the Monte Carlo noise is suppressed 
by a high number of scenarios, the n max = 3 case exhibits better performance in terms of accuracy. The 
agreement between the two methods in case n max = 3 and N = 10 8 is remarkably good. 

For both n max = 2 and n ma x = 3 cases the deviations from the a ~ Monte Carlo error convergence 

rule have been studied. Figure [2] shows the dependency of the a ■ on the number of simulations. 

For the n max = 3 case the product a ■ \f~N remains constant up to N = 10 8 , which implies that the 
difference between the Monte Carlo and analytical estimates is mainly due to the noise of the Monte 
Carlo simulation. Thus, the accuracy of the analytical approximation in case n max = 3 exceeds the 
accuracy of Monte Carlo simulation using 10 s scenarios. 

Computer time spent was 13 seconds for n max = 3 analytical and 16 hours for N = 10 s Monte Carlo 
calculations. Overall, the proposed analytical approximation to variance-covariance based risk allocation 
demonstrated exceptionally good performance both in terms of accuracy and computer time. 

3.4 Accuracy of the approximation 

The high accuracy of the proposed approximation has been demonstrated by benchmarking the analytical 
estimates against Monte Carlo simulation. However, some explanations are necessary to address possible 
scepticism regarding the convergence properties of the expansion (|2.11|) . 

Indeed, approximating a pairwise correlation by the first three terms of the above mentioned series 
expansion may lead to a rather rough estimate. Asset correlation pij of two borrowers with high systematic 
components of asset returns (rj in (|2.5p ) from the same country/industry cluster (i.e. the same factor 
loadings {Pi}) may be higher than 0.5, which results in a poor convergence of (|2.1ip . However, two 
important mitigating factors contribute to the overall accuracy of the approximation. 

First, according to (|2.3|) . the risk contribution of is proportional to a sum of covariances of the facility 
i with all other facilities in a portfolio. The accuracy in this case depends on the average pairwise asset 
correlations which are quite smaller (compared to 0.5) in realistic portfolios due to country/industry 
diversification and existence of facilities with low systematic components of asset returns. 

Second, contributions to the higher order (n > 3) parameters k in (|2.15p have different 
signs depending on the valuation function Uj(e). This pd diversification effect is obvious in a default-only 
case (|2.21j) and is caused by oscillating nature of higher order Hermite polynomials. As a consequence, 
the higher order contributions decrease and overall convergence properties of (|2.17p arc improved. 

Both of the above mentioned factors contributing to the accuracy of the proposed approximation are 
significant if the portfolio consists of large number of loans with different characteristics (country/industry 
association, probability of default, size of systematic part of asset return, etc). In other words, the more 
realistic the portfolio, the better the accuracy of the approximation. 

a is a standard deviation of the relative difference between analytical and simulation-based estimates of risk contributions 
<r c across the portfolio. 
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4 Summary 



Variance-covariance based risk allocation method has been proposed. The method was applied to a " real 
life" credit portfolio model and tested on a "real life" portfolio. It was shown that the proposed analytical 
approximation is superior both in terms of speed and accuracy to the traditional Monte Carlo simulations. 

The proposed routine may be especially appealing for the purpose of risk adjusted loan pricing since 
it allows to compute accurate and statistical noise-free estimates of capital charges. Once the portfolio- 
specific parameters (|2.15[) are calculated, the risk contribution <j c of a loan can be computed using (|3.10|) 
and the capital charge can be assigned according to (|2.4| . 
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